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ABSTRACT 

A pulsar timing system has been operating in the 18-cm band at the Urumqi As- 
tronomical Observatory 25-m telescope since mid-1999. Frequency resolution allowing 
dedispersion of the pulsar signals is provided bya2xl28x2.5 MHz filterbank/digitiser 
system. Observations of 74 pulsars over more than 12 months have resulted in updated 
pulsar periods and period derivatives, as well as improved positions. Comparison with 
previous measurements showed that the changes in period and period derivative tend 
to have the same sign and to be correlated in amplitude. A model based on unseen 
glitches gives a good explanation of the observed changes, suggesting that long-term 
fluctuations in period and period derivatives are dominated by glitches. In 2000 July, 
we detected a glitch of relative amplitude Av/v ~ 24 x 10~^ in the Crab pulsar. The 
post-glitch decay appears similar to other large Crab glitches. 
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1 INTRODUCTION 

Pulsars are remarkable celestial clocks. They act as probes 
for many difFerent studies, including investigations of the in- 
terstellar medium, stellar and binary system evolution, rel- 
ativistic astrophysics and condensed matter. Identified as 
rotating neutron stars, they have an observed pulse period 
range of 0.0016 to 8.5 s. Because of the huge moment of in- 
ertia (~ 10*^ g cm^) and compact nature of a neutron star, 
the rotation periods are very stable. However, they are not 
completely constant, increasing gradually as the pulsar loses 
energy in the form of relativistic particles and electromag- 
netic radiation with small variations due to changes in the 
stellar structure. Pulse timing observations can be used to 
accurately measure the pulse periods and these variations, 
including effects related to the orbital motion for binary pul- 
sars, precise positions and proper motions and even annual 
parallax for some pulsars. 

If we express the pulsar rotation frequency as a Taylor 
series, the pulse phase at time t is given by 



(1) 



where the rotation frequency v = 1/P, where P is the pulse 
period, i/ and i> are the first and second time derivatives of 
the pulse frequency, and t is the time in a frame which is 
inertial with respect to the pulsar, normally taken to be the 
solar-system barycenter frame. The differences between ob- 
served pulse arrival times (TO As) and those predicted from 
Equation |l|, are known as timing residuals. Timing residu- 
als will deviate from zero if the model is not accurate, for 
example, due to errors in the frequency parameters or irregu- 
larities in the pulse period. Given a set of TO As from timing 
observations at different epochs t, corrections to the rotation 
and astrometric parameters can be obtained by performing 
a least-square fit to minimise the residuals. Improved po- 
sitions, with an accuracy of the order 0.1 arcsec, can be 
obtained from TOAs spanning a year or more, because an 
error in position gives rise to an annual term in the timing 
residuals due to the orbital motion of the Earth. 

Two kinds of timing irregularities are observed in pul- 
sars. The first is timing noise, characterised by random pe- 
riod changes with the fractional amplitude about 10~^ s or 
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less with time scales of days, months or years ( Cordes & 
Downs 1985; Arzoumanian et al. 1994). The second kind of 
timing irregularity is known as a 'glitch', where the pulse 
frequency suddenly increases with the fractional amplitude 
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^vjv t ypically in the range 

10-9 „ 10 -6 ( |Lyne, Shemar fc| 
Grahan ji- Smith 200(]| ; [Wang et al. 2000[ ). Glitches are unpre- 
dictable but typically occur at intervals of a few years in 
young pulsars. Coincident with the glitch, there is often an 
increase in the magnitude of the frequency derivative, typi- 
cally by ~ 1 per cent, which sometimes has an approximately 
exponential decay with a time scale of weeks to years. 

Because of these unpredictable changes, period param- 
eters lose their validity with time. However, for some pul- 



sars, p arameters have not bee n updated for 30 years (Tay- 
lor, Manchester fc Lyne 1993), and hence it is desirable to 



improve them. Moreover, to study the properties of timing 
noise and glitches, frequent timing observations are required. 

The 25-m radio telescope at Nanshan, operated by 
Urumqi Astronomical Observatory (UAO), is ideal for such 
applications. It is mainly used as a VLBI station and has 
receivers for six wavelength bands: 1.3 cm, 3.6/13 cm, 6 cm, 
18 cm and 92 cm. The telescope is available for pulsar tim- 
ing for about two days every week, and so frequent obser- 
vations are possible. To overcome the undesirable effects 
of high Galactic background emission and of propagation 
in the interstellar medium, especially scattering at longer 
wavelength, we decided to use the 18-cm band for the pul- 
sar timing observations. A timing system based on a room- 
temperature receiver has been operational since mid-1999. 
In Section ^ of this paper, the receiver and data acquisition 
system are introduced. Observations and data analysis are 
discussed in Section ^. In Section ^, we present new mea- 
surements of positions, discuss period and period derivative 
changes and introduce a glitch detected in the Crab pulsar. 
We discuss and summarise our work in Section |^. 



2 RECEIVER AND DATA ACQUISITION 
SYSTEM 

A pulsar timing system has been established at the Nanshan 
Station of Urumqi Astronomical Observatory, located near 
the geographic center of Asia, with longitude 87° and lati- 
tude -1-43° . The timing system operates in the 18-cm band 
and is outlined in Fig. hi For this band, the telescope has 
cassegrain optics and uses a horn feed receiving orthogonal 
circular polarisations. The receiver has dual-channel, room- 
temperature pre-amplifiers with center radio frequency (RF) 
of 1540 MHz and total bandwidth of 320 MHz. The system 
temperature is approximately 100 K. The polarisations are 
split by an ortho-mode transducer (OMT) at the end of the 
feed, amplified and then down-converted to intermediate fre- 
quency (IF) in the range 80-400 MHz using a local oscillator 
(LO) at 1300 MHz. After conversion, the signals are fed to a 
filterbank system which has 128 2.5-MHz channels for each 
polarisation. 

Each of the 256 analogue signals from the filterbank 
pass through a 0.2 Hz high-pass filter and are then integrated 
and 1-bit digitised. Digitised data are grouped in blocks of 
32 bits for transfer to the on-line computer. The digitiser 
uses 1-sec and 5-sec pulses derived from the Observatory H- 
maser to start the data conversion at a known time. The 
observation is started by the 1-sec pulse and the start time 
is checked using the 5-sec pulse. The sampling interval is 
adjustable as a multiple of 1 /is, clocked by a 1 MHz reference 
signal, also derived from the Observatory H-maser. The first 
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Figure 1. Diagram of the pulsar timing system at Urumqi Astro- 
nomical observatory. The center frequency is 1540 MHz and the 
received bandwidth and the received band is defined by bandpass 
filters (BPF) after the RF amplifiers. The ortho-mode transducer 
(OMT) and the receiver preamplifiers are currently at room tem- 
perature but will soon be replaced by a cryogenically cooled sys- 
tem. De-dispersion is provided by a 2 X 128 X 2.5 MHz filterbank 
and digitiser. The digitiser is controlled by a PC computer using 
National Instruments interface cards. 



word of each sample is a counter which is checked by the data 
acquisition program. A GPS time-transfer system is used to 
align the Observatory clock with UTC. 

Signals from each channel are recorded by a data ac- 
quisition system based on a PC operating under Windows 
NT. The online program (Pulsar Timing Data Acquisition) 
is written in Visual C-I--I-. It has threads to perform sev- 
eral tasks at the same time, including setting the digitiser 
parameters and data sampling, unpacking data and check- 
ing time synchronisation, folding of channel data and saving 
data to disk, dedispersing and displaying of folded profiles 
and communicating with the computer which controls the 
telescope. The program reads in the date and UTC from 
the clock. A run-line is set high before the chosen start time 
and the sampling is triggered by the rising edge of the next 
1-sec pulse. Unpacked data for each channel are folded at 
the apparent pulsar period, updated at 120-s intervals, to 
form sub-integrations of typical duration 240 s. Predicted 
topocentric periods are calculated using 'polyco' files pro- 
duced by the pulsar timing program TEMPOQ Folded data 



t See http://www.atnf.csiro.au/rosearch/pulsar/timing/tempo/ 
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Figure 2. Pulse profiles of PSR B1933+16. The upper part shows 
the profile in each of the 128 2.5MHz channels, the lower part is 
the sum of all the channels after appropriate shifts for dispersion 
have been applied. 



are saved to disk at every sub-integration time. The sequence 
counter and synchronisation flag produced by the digitiser 
are checked every sample and on each 5 second respectively. 
Folded profiles, including each channel, the sum of each po- 
larisation and total power, are displayed for monitor pur- 
poses (See Fig. gI) and are refreshed about once per second. 



3 OBSERVATIONS AND DATA ANALYSIS 

The sensitivity of pulsar observations is affected by re- 
ceiver, spillover and Galactic background noise (Tree, Tspi 
and Tbkg), antenna gain (G), receiver bandwidth (Ai^), in- 
tegration time (tint), effective pulse width (W) and pulse 
period (P). The limiting mean flux density is given by: 



Q/3(Tbkg + Jspl + Tree) 

iiit 



W 



P-W 



(2) 



where a ~ 5 is the threshold signal-to-noise ratio for timing, 
P ~ 1.5 is digitisation and other processing losses, Tbkg -|- 
Tspi + Tree ~ 100 K, G ~ 0.1 K Jy"\ iV^ = 2 is the number 



The pulsar catalogue ( Taylor, Manchester & Lyne 199J ) 
gives flux densities, S, for most pulsars at 400 MHz and 
1400 MHz. We first selected pulsars that are north of decli- 
nation —40° and have 51400 > 5min. For those where there 
are no 1400 MHz measurements, we assume a power law 
spectrum S oc where / is the observing frequency and 
Q — —1.7, to estimate 51400 from S'4oo- AH these pulsars 
were observed from October 1999 to December 1999. From 
December, only those pulsars where significant emission was 
detected in 16 minutes were observed. There were 74 such 
pulsars and these generally had ^1400 > 4 mjy. All of these 
pulsars are 'normal' (non-millisecond) pulsars. Their posi- 
tions on a P — P diagram are shown in Fig. ^ Observations 
are made once or twice per week with most being of 16 
minutes duration. In this paper we present observations for 
these 74 pulsars up to February 2001. 

After each observation, the data are moved to a Linux 
system and compressed using the program TREDUCE (sup- 
ported by Swinburne University of Technology and ATNF) 
to form the archive files. These archive flies consist 8 sub- 
bands in each polarisation and 4 sub-integrations. The 
archives are first dedispersed relative to the center frequency 
and summed in time by TREDUCE to form one mean pulse 
profile per observation. TOAs are obtained using an itera- 
tive process. First a profile with good signal-to-noise ratio is 
adopted as a 'standard profile' or template. This template 
is cross-correlated with the observed mean pulse profiles to 
produce TOAs. TEMPO is used to reduce TOAs to arrival 
times at inflnite frequency at the Solar syste m barycentre us - 
ing the JPL Solar system ephemeris DE200 (Standish 1982), 



of polarisations and Au — 320 MHz ( Dewey et al. 1985 ). For 
a typical pulsar with a duty cycle W/P ~ 0.05, Equation ^ 
gives a sensitivity of Smin ~ 2.2 mJy for tint ~ 16 mins. 



compute timing residuals and do a least-square flt to de- 
termine a preliminary set of improved model parameters. 
Using these parameters, all observed proflles are summed 
to produce a high signal-to-noise ratio mean pulse profile 
which is then adopted as a standard profile for subsequent 
TEMPO analyses. As examples. Fig. U shows six final stan- 
dard profiles. For the mode-changing pulsar PSR B0329-(-54, 
two standard profiles are formed for the different modes and 
TOAs are obtained using the appropriate profile. 



4 RESULTS 

In this section we present the results of timing analyses of 
74 pulsars. These results are compared with previously pub- 
lished data from the pulsar catalogue. In addition we checked 
recent publications on pulsar timing and astrometry to ob- 
tain the best periods, period derivatives and positions. We fit 
for position, period and period first derivative for all pulsars 
except as noted below. The epochs of periods and positions 
are fixed at MJD 51700, near the center of our data span. 
PSRs B1754-24 and B1800-21 lie near the ecliptic plane 
and so their declinations cannot be determined from timing 
fits. The Crab pulsar suffered a glitch during the year and 
is discussed separately. 



4.1 Positions 

Pulsar positions from the catalogue and this work are 
present in Table |^; superscripts 'c' and 'o' stand for cat- 
alogue value and observed value respectively. Parameter un- 
certainties in this and subsequent tables are given in paren- 
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Figure S. P — P diagram of known pulsars. Pulsars observed 
at UAO are marked as stars. Binary pulsars are indicated by 
a circle around the point and lines of constant surface dipole 
magnetic field and characteristic age are shown. The spin-up line 
for millisecond pulsars is indicated by a dashed line. 
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Figure 4. Examples of standard pulse profiles obtained from the 
observations which are used as templates to form the TOAs. The 
full pulse period is shown for each profile. 



theses and refer to the last quoted digit. For the observed 
parameters, they are the 1-a errors given by TEMPO. Ref- 
erences for catalogue positions are given in column 5 and the 
data span of our observations is given in the last column. 

Among these pulsars, PSR B1754— 24 had a very large 
sinusoid timing residual before fitting for position. Since it 
lies close to the ecliptic plane, timing measurements can- 
not determine an accurate declination. The best previously 
published measurements are R.A.(J2 000) = 17''57™40''(60) 
( [Taylor, Manchester fc Lyne 199^) and Dec.(J200q ) = 



24°21'57"(30) (Vivekanand, Mohanty fc Saher 1983). In 



our analysis, we fix the declination at the catalogue value 
and fit for R.A. along with the other timing parameters. The 
error of R.A. was estimated by taking the quadrature sum of 
the TEMPO R.A. error at the nominal declination and half 
of the difference in the R.A. in TEMPO fits with the decli- 
nation held at each end of its uncertainty range. This gives 
the greatly improved value R.A.(J2000) = 17''57™29':32(3). 
Similarly, for PSR B1800-21 we obtained R.A.(J2000) = 
18'^03™51'!15(3), significantly different from the catalogue 
value. 

In principle, the parameters given in Table ^ can be 
used to estimate the pulsar proper motions and (two- 
dimensional) space velocities. However, with only one year of 
data, it is impossible to separate the effects of proper motion 
and long-period timing noise, leading to possible systematic 



errors in the derived proper motions. We therefore defer this 
analysis to a time when more data are available. 



4.2 Period Parameters 

Improved periods and period derivatives at epoch MJD 
51700 are given in Table ^ for all observed pulsars except 
the Crab pulsar. Column 5 gives the number of TOAs and 
column 6 the timing residual in microseconds. 

By comparing the observed and catalogue values, we 
obtain the period and period derivative changes, AP and 
AP given in columns 7 and 8 respectively. The time inter- 
val from the epoch of the catalogue period to MJD 51700 
is given in column 9. There are seven pulsars for which 
glitches have been detected in previous observations: PSRs 
B0355+54, B0525+21, B0531+21, B1508+55, B1737-30, 



B1800-21 and B1830- 08 (Lyne, Shemar fc Graham-Smith 



200C ; Wang et al. 2000). We quote the latest published post- 



glitch rotation parameters for comparison. The Crab pulsar 
(PSR B0531-I-21) glitched in July 2000; results for it are 



discussed in Section 4.3 



Significant differences between the observed and pre- 
dicted periods are observed for most pulsars. Some pul- 
sars, such as PSRs B0628-28, B0844-35, B1737-39, 
B1834— 10 and B1845— 01, the period changes are posi- 
tive and larger than 10~* s. In contrast, PSRs B1700— 32, 
B1706-16, B1754-24 and B1822-09 have periods about 
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Table 1. Pulsars positions from catalogue and this work. 
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-40: 


:53: 
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10 


17: 


:07: 


:21, 


,728(12) 


-40: 


:53: 
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51573- 


-51925 


1709-1640 
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:40: 


:58.4(3) 
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:45.6(4) 


51506- 


-51967 
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:40: 
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.7(1) 


-30: 
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11 


17: 


:40: 


:33, 


,753(16) 


-30: 


:15: 
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51549- 
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1741-3927 
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17: 
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:18, 


.04(3) 


-39: 
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17: 


:41: 
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:27: 


:38.41(20) 


51549- 
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17: 
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:29, 
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:51, 
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11 


18: 


:03: 


:51, 
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:38, 
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:23, 
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:20: 


:52, 
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:27: 
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1 


18: 


:20: 


:52, 


,620(5) 


-04: 


:27: 


:37.1(4) 
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-51967 


1824-1945 


1821-19 


18: 


:24: 


:00, 


.45(4) 


-19: 


:45: 


:51(8) 


8 


18: 


:24: 


:00, 


,449(6) 


-19: 


:45: 


:49.9(12) 


51512- 


-51966 


1825-0935 


1822-09 


18: 


:25: 


:30, 


.596(6) 


-09: 


:35: 


:22.8(4) 


2 


18: 


:25: 


:30, 


,599(7) 


-09: 


:35: 


:21.9(6) 


51506- 


-51967 


1829-1751 


1826-17 


18: 


:29: 


:43, 


.12(1) 


-17: 


:51: 


:03(2) 


7 


18: 


:29: 


:43, 


,130(10) 
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;40, 
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-08: 


:27: 
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14 


18: 


:33: 


:40, 
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-08: 


:27: 


;31. 50(15) 


51597- 


-51960 


1835-1106 




18: 


:35: 


:18, 


.287(2) 


-11: 


:06: 


:15.1(2) 


15 


18: 


:35: 


:18, 


,17(3) 


-11: 


:06: 


:19.3(13) 
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-51966 


1836-1008 


1834-10 


18: 


:36: 


:53, 
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-10: 
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:09(5) 


16 


18: 


:36: 


:53, 


,911(5) 


-10: 


:08: 


:09.8(4) 
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-51966 


1840+5640 


1839+56 
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:57.9(3) 


51500- 
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" 1. Downs & Reichley, 1983 2. Arzoumanian et al., 1994 3. Fomalont et al., 1992 4. Fomalont et al., 1997 5. Dewey et al., 
1988 6. McNamara, 1973, 7. Siegman, Manchester & Durdin, 1993 8. Newton, Manchester & Cooke, 1981 9. Manchester 
& Taylor, 1981 10. Johnston et al., 1995 11. Braun, Goss and Lyne, 1989 12. Taylor, Manchester and Lyne, 1993 13. 
Vivekanand, Mohanty & Salter, 1983 14. Clifton et al., 1992 15. Manchester et al., 1996 16. Backus, Taylor & Damashek, 
1982 17. GuUahorn & Rankin, 1978 18. Foster, Backer & Wolszczan, 1990 
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Table 1. -Continued. 



PSR J 


PSR B 






Ref. 






Data Span 






fh m s \ 
V / 






/h m s "1 
V ) 


\ } 


MJD 


1900-2600 


1857-26 


19:00:47.60(1) 


-26:00:43.1(3) 


3 


19:00:47.561(18) 


-26:00:39(3) 


51573- 


-51967 


1913-0440 


1911-04 


19:13:54.18(1) 


-04:40:47.6(4) 


9 


19:13:54.184(5) 


-04:40:47.6(3) 


51550- 


-51966 


1917+1353 


1915+13 


19:17:39.784(2) 


+13:53:57.06(9) 


2 


19:17:39.790(4) 


+13:53:56.80(9) 


51573- 


-51966 


1921+2153 


1919+21 


19:21:44.798(3) 


+21:53:01.83(8) 


9 


19:21:44.80(3) 


+21:53:02.8(5) 


51600- 


-51966 


1932+1059 


1929+10 


19:32:13.900(2) 


+10:59:31.99(7) 


2 


19:32:13.9430(15) 


+10:59:32.52(6) 


51550- 


-51960 


1935+1616 


1933+16 


19:35:47.835(1) 


+16:16:40.59(2) 


1 


19:35:47.8249(11) 


+16:16:40.03(4) 


51550- 


-51966 


1946+1805 


1944+17 


19:46:53.043(4) 


+18:05:41.59(9) 


17 


19:46:53.05(3) 


+18:05:41.5(8) 


51560- 


-51966 


1948+3540 


1946+35 


19:48:25.037(2) 


+35:40:11.28(2) 


17 


19:48:25.006(3) 


+35:40:11.03(6) 


51500- 


-51966 


1955+5059 


1953+50 


19:55:18.90(6) 


+50:59:54.2(6) 


16 


19:55:18.736(6) 


+50:59:55.76(10) 


51513- 


-51966 


2002+4050 


2000+40 


20:02:44.04(3) 


+40:50:54.7(3) 


5 


20:02:44.047(10) 


+40:50:53.99(19) 


51560- 


-51966 


2013+3845 


2011+38 


20:13:10.49(3) 


+38:45:44.8(3) 


5 


20:13:10.384(7) 


+38:45:43.08(16) 


51549- 


-51966 


2018+2839 


2016+28 


20:18:03.851(2) 


+28:39:54.26(3) 


1 


20:18:03.832(4) 


+28:39:54.26(9) 


51549- 


-51966 


2022+2854 


2020+28 


20:22:37.079(3) 


+28:54:23.45(3) 


18 


20:22:37.078(3) 


+28:54:22.91(6) 


51547- 


-51967 


2022+5154 


2021+51 


20:22:49.900(2) 


+51:54:50.06(2) 


1 


20:22:49.843(10) 


+51:54:50.46(16) 


51547- 


-51968 


2048-1616 


2045-16 


20:48:35.472(4) 


-16:16:44.45(8) 


3 


20:48:35.46(8) 


-16:16:38(4) 


51547- 


-51960 


2108+4441 


2106+44 


21:08:20.48(1) 


+44:41:48.8(1) 


2 


21:08:20.521(8) 


+44:41:49.37(16) 


51547- 


-51966 


2113+4644 


2111+46 


21:13:24.29(1) 


+46:44:08.7(1) 


2 


21:13:24.377(10) 


+46:44:09.1(3) 


51547- 


-51966 


2157+4017 


2154+40 


21:57:01.82(1) 


+40:17:45.9(1) 


2 


21:57:01.879(13) 


+40:17:45.8(4) 


51547- 


-51966 


2219+4754 


2217+47 


22:19:48.136(4) 


+47:54:53.83(4) 


2 


22:19:48.115(10) 


+47:54:53.73(12) 


51506- 


-51966 


2257+5909 


2255+58 


22:57:57.711(4) 


+59:09:14.95(3) 


2 


22:57:57.756(9) 


+59:09:14.93(9) 


51547- 


-51966 


2313+4253 


2310+42 


23:13:08.571(6) 


+42:53:12.98(3) 


3 


23:13:08.622(10) 


+42:53:13.12(19) 


51547- 


-51967 


2321+6024 


2319+60 


23:21:55.19(4) 


+60:24:30.7(3) 


2 


23:21:55.16(3) 


+60:24:30.58(16) 


51499- 


-51966 


2326+6113 


2324+60 


23:26:58.704(5) 


+61:13:36.50(3) 


2 


23:26:58.693(8) 


+61:13:36.29(8) 


51560- 


-51966 


2354+6155 


2351+61 


23:54:04.71(2) 


+61:55:46.8(1) 


2 


23:54:04.743(16) 


+61:55:46.67(12) 


51547- 


-51966 



" 1. Downs & Reichley, 1983 2. Arzoumanian et al., 1994 3. Fomalont et al., 1992 4. Fomalont et al., 1997 5. Dewey et al., 
1988 6. McNamara, 1973, 7. Siegman, Manchester & Durdin, 1993 8. Newton, Manchester & Cooke, 1981 9. Manchester 
& Taylor, 1981 10. Johnston et al., 1995 11. Braun, Goss and Lyne, 1989 12. Taylor, Manchester and Lyne, 1993 13. 
Vivekanand, Mohanty & Salter, 1983 14. Clifton et al., 1992 15. Manchester et al., 1996 16. Backus, Taylor & Damashek, 
1982 17. GuUahorn & Rankin, 1978 18. Foster, Backer & Wolszczan, 1990 



10 * s less than predicted. The observed period change in 
PSR B1737-30 is huge, with AP ~ 1.2 x 10"'' s, much 
larger than for other pulsars. This pulsar giitched frequently 
during MJD 46301-49 400, with nine glitches being detected 
( Shemar fc Lyne 1996 ), i.e., on average, one glitch every 11 



months It has been more than six years since Shemar & 



Lyne's work until our observations started, so it is probable 
that there were several glitches during this interval. 

Fig. ^ shows a histogram of the observed period changes 
excepting PSR B1737— 30 and the Crab pulsar. Also, pul- 
sars are not included if the uncertainties in AP or AP 
are larger than 20 ns or 30 x 10~^* respectively, eliminat- 
ing PSRs B0525+21, B0844-35, B1700-32, B1754-24 and 
B1800— 21. With these exceptions, most period changes are 
small, with |AP| ;^3 ns. 

We find an interesting correlation between AP and AP, 
illustrated in Fig. ^. This plot shows a strong tendency for 
AP and AP to have the same sign and to be correlated in 
amplitude. For random period noise we would expect these 
two quantities to be uncorrelated. However, glitches involve 
systematic changes in both P and P and we investigate this 
further. 

In most cases, the pulse frequency after a glitch is well 
modelled by the relation 



where vo{t) is the value of v extrapolated from before the 
glitch, /S.Vg is the total frequency change at the time of the 
glitch, Q is the fractional part of A/Zg which decays expo- 
nentially, Td is the decay time constant and A;> p is the per- 



mane nt change in v at the time of the glitch (Wang et al 



200C). The increment in frequency derivative is given by: 



AO{t) = Ji^t!l. exp(-t/rd) + A!>p. 



Td 



(4) 



v{t) = vo{t) + Ai/g[l - Q(l - exp(-t/rd)] + Ai-pt, 



(3) 



We use these relations to compute P and P at an ar- 
bitrary time t in the past, and compare them with the val- 
ues at the current time to, assuming a glitch of magnitude 
Afg = lO"*' at time t^-t = 400 days. Values of AP = Pq-P 
and AP = Pq — P are shown as a function of to — i in Fig. ^ 
If a pulsar giitched before the previous observation, both 
AP and AP will be negative; while if the glitch happened 
between the observations, we have AP < and AP < 
for negative APp = -Ai/p/i/^, and AP > and AP > for 
positive APp and large fractional decay Q. These predictions 
are consistent with the observations shown in Fig. ^ The 
last case represents the observed sit uation in the Crab pul- 
sar ( Lyne, Pritchard fc Smith 1993| ). For positive APp and 
small Q, we can have AP < and AP > 0. Furthermore, in 
this model there is no case with AP > and AP < 0, consis- 
tent with the observed deficit in the lower-right of Fig. ^ For 
Q > in a glitch before the first observation, the changes in 
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Table 2. Pulsar period parameters at MJD 51700. 



PSR J 



PSR B 



P 

(s) 



P 

(10-15) 



Res. 



AP 

(10-9 s) 



AP 

(10-18) 



AT Ref.[ 

(yr) 



0034—0721 


0031- 


-07 


0.94295117325(7) 


0.395(13) 


25 


854 


—0.03(87) 


— 13(13) 


30.1 


1 


0139+5814 


01364 


-57 


0.272452862190(3) 


10.7040(6) 


60 


268 


3.07(9) 


3.4(6) 


9.1 


2 


0141+6009 


0138+59 


1.22294860175(7) 


0.417(17) 


41 


756 


0.4(7) 


26(17) 


27.1 


3 


0332+5434 


0329+54 


0.714520624711(5) 


2.0489(10) 


102 


469 


— 1.19(8) 


-0.7(10) 


30.3 


1 


0358+5413 


03554 


-54 


0.1563832093117(10) 


4.3954(3) 


61 


342 


0.66(10) 


5.0(5) 


14.2 


4 


0452—1759 


0450- 


-18 


0.548940421758(19) 


5.753(3) 


54 


546 


— 1.31(18) 


-4(3) 


13.4 


5 


0454+5543 


0450+55 


0.340729803314(5) 


2.3742(8) 


59 


321 


1.95(6) 


8.6(8) 


9.0 


2 


0528+2200 


0525+21 


3.74553060709(15) 


40.01(3) 


52 


893 


2(11) 


— 16(36) 


26.4 


4 


0543+2329 


05404 


-23 


0.245978510226(3) 


15.4205(7) 


65 


568 


—0.764(18) 


-3.8(7) 


9.1 


2 


0612+3721 


06094 


-37 


0.297982336991(18) 


0.059(3) 


21 


293 


0.23(6) 


0.05(300) 


15.3 


6 


0630—2834 


0628- 


-28 


1.24442173630(8) 


7.195(14) 


53 


1171 


11(6) 


88(16) 


20.7 


7 


0738—4042 


0736- 


-40 


0.374919985032(6) 


1.6161(9) 


44 


171 


-1.3(3) 


2.0(9) 


25.0 


7 


0742—2822 


0740- 


-28 


0.166765742577(5) 


16.8177(8) 


62 


575 


3.20(6) 


5.2(8) 


9.1 


2 


0814+7429 


08094 


-74 


1.29224148381(11) 


0.153(16) 


57 


1119 


0.03(26) 


— 15(16) 


9.1 


2 


0820-1350 


0818- 


-13 


1.23813005246(4) 


2.097(6) 


58 


644 


-0.24(19) 


-8(6) 


29.3 


1 


0826+2637 


08234 


-26 


0.530661285865(12) 


1.6978(19) 


66 


588 


— 1.61(12) 


— 11.6(20) 


9.1 


2 


0837—4135 


0835- 


-41 


0.751623617646(8) 


3.5393(12) 


33 


165 


-1.8(3) 


—7.6(13) 


13.4 


5 


0846—3533 


0844- 


-35 


1.1160975763(3) 


1.62(5) 


30 


1046 




48(53) 


22.3 


8 


0922+0638 


09194 


-06 


0.430623569113(12) 


13.6924(19) 


54 


521 


— 1.18(15) 


—28.1(20) 


9.0 


2 


0953+0755 


09504 


-08 


0.253065270746(6) 


0.2288(10) 


69 


540 


0.62(10) 


-0.4(10) 


27.9 


9 


1136+1551 


11334 


-16 


1.18791477322(3) 


3.733(5) 


73 


661 


0.60(7) 


0.6(48) 


27.5 


9 


1239+2453 


12374 


-25 


1.38244953180(4) 


0.967(7) 


67 


677 


—0.15(15) 


6(7) 


9.1 


2 


1509+5531 


15084 


-55 


0.739682698169(16) 


4.996(3) 


58 


512 


—2.68(9) 


-12(3) 


9.1 


2 


1645—0317 


1642- 


-03 


0.387690495052(10) 


1.7924(15) 


67 


595 


— 1.04(10) 


11.4(15) 


30.3 


1 


1703-3241 


1700- 


-32 


1.21178519124(4) 


0.660(7) 


57 


709 


-36(31) 


-130(97) 


10.1 


2 


1705—1906 


1702- 


-19 


0.298988488005(8) 


4.1370(14) 


55 


603 


—0.244(20) 


-1.3(14) 


9.2 


2 


1707—4053 


1703- 


-40 


0.581016707860(18) 


1.925(4) 


23 


280 


0.19(7) 


5(4) 


9.1 


2 


1709—1640 


1706- 


-16 


0.65305652375(3) 


6.240(4) 


50 


675 


— 11.99(12) 


—69(4) 


9.1 


2 


1721—3532 


1718- 


-35 


0.280424567967(9) 


25.1816(14) 


52 


497 


—4.28(9) 


— 14.5(15) 


9.1 


2 


1722-3207 


1718- 


-32 


0.477157619208(5) 


0.6475(9) 


57 


378 


-0.25(4) 


0.4(9) 


9.1 


2 


1740—3015 


1737- 


-30 


0.60676381008(3) 


466.312(5) 


48 


893 


— 1181(9) 


291(54) 


6.7 


4 


1741—3927 


1737- 


-39 


0.512211401067(8) 


1.9313(14) 


44 


319 


68(3) 


122(5) 


22.3 


8 


1745—3040 


1742- 


-30 


0.367430422265(6) 


10.6696(10) 


54 


327 


0.86(3) 


4.5(10) 


9.1 


2 


1752—2806 


1749- 


-28 


0.562561300009(9) 


8.1271(16) 


62 


550 


2.0(7) 


5(3) 


9.5 


5 


1757-2421 


1754- 


-24 


0.234102570587(4) 


12.9146(5) 


51 


373 


-54(192) 


-100(300) 


20.2 


10 


1803—2137 


1800- 


-21 


0.133641979388(15) 


134.150(3) 


36 


2378 


5(6) 


206(200) 


1.9 


11 


1807—0847 


1804- 


-08 


0.163727380271(3) 


0.0287(5) 


58 


725 


0.09(7) 


0.03(48) 


21.4 


3 


1818—1422 


1815- 


-14 


0.291489038008(5) 


2.0399(9) 


49 


316 


0.9(3) 


3.0(11) 


12.9 


12 


1820—0427 


1818- 


-04 


0.598078699467(10) 


6.3379(20) 


60 


744 


-6.1(5) 


0.2(21) 


30.3 


1 


io24 — iy40 


1821- 


-19 


0. 189335814049(3 j 


0.2394(0 ) 


63 


519 


2.3(6) 


14.0(14 ) 


13.4 


5 


1825-0935 


1822- 


-09 


0.768994380569(16) 


52.336(4) 


59 


909 


-33.5(3) 


-29(4) 


9.1 


2 


1829-1751 


1826- 


-17 


0.307133849665(11) 


5.5414(16) 


50 


420 


-4.17(13) 


-20.6(16) 


13.4 


5 


1832-0827 


1829- 


-08 


0.647305181037(9) 


63.8874(20) 


55 


713 


-3.70(17) 


-5.8(20) 


12.9 


12 


1833-0827 


1830- 


-08 


0.0852852156141(15) 


9.1721(3) 


39 


145 


0.30(5) 


2.6(3) 


10.0 


4 


1835-1106 






0.165910996503(20) 


20.656(3) 


27 


915 


3.896(20) 


44(3) 


6.2 


13 


1836-1008 


1834- 


-10 


0.562713612837(7) 


11.8013(14) 


58 


618 


15(5) 


26(7) 


21.4 


3 


1840+5640 


1839+56 


1.65286223816(14) 


1.51(3) 


46 


1812 


-0.01(88) 


10(24) 


9.1 


2 


1847-0402 


1844- 


-04 


0.59778239230(3) 


51.704(4) 


50 


708 


-2.3(5) 


-10(4) 


13.4 


5 


1848-0123 


1845- 


-01 


0.659432812074(9) 


5.2479(18) 


58 


709 


11.6(7) 


29(3) 


13.4 


5 


1900-2600 


1857- 


-26 


0.61220925419(4) 


0.193(6) 


49 


788 


0.20(6) 


-11(6) 


9.1 


2 



° 1. Manchester & Taylor, 1981 2. Arzoumanian et al., 1994 3. Backus, Taylor & Damashek, 1982 4. Shemar 
& Lyne, 1996 5. Siegman, Manchester & Durdin, 1993 6. Dewey et al., 1988 7. Manchester et al. 1983 8. 
Newton, Manchester & Cooke, 1981 9. GuUahorn & Rankin, 1978 10. Vivekanand, Mohanty & Salter, 1983 
11. Wang et al., 2000 12. Clifton et al., 1992 13. Manchester et al., 1996 14. Foster, Backer & Wolszczan, 1990 
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Table 2. Pulsar period parameters at MJD 51700. 



PSR J 


PSR B 


P 


P 




Res. 


AP 


AP 


AT 


Ref.F| 
U 








(s) 


(10-15) 




(/js) 


(10-9 s) 


(10-18) 


(yr) 




1913-0440 


1911- 


-04 


0.825937582927(13) 


4.055(3) 


47 


301 


-1.37(20) 


-14(3) 


30.3 


1 


1917+1353 


1915- 


hl3 


0.1946321872566(13) 


7.19502(18) 


49 


256 


-0.150(18) 


-3.22(19) 


9.1 


2 


1921+2153 


1919- 


1-21 


1.33730247470(8) 


1.358(8) 


26 


1151 


-0.02(9) 


10(8) 


30.1 


1 


1932+1059 


1929- 


hlO 


0.2265181529285(15) 


1.1618(3) 


39 


120 


0.39(4) 


5.2(3) 


9.1 


2 


1935+1616 


1933- 


hl6 


0.3587411416716(16) 


6.0023(3) 


42 


84 


-0.91(5) 


-1 4(3) 


25.8 


9 


1946+1805 


1944- 


hl7 


0.44061848299(3) 


0.020(3) 


54 


2234 


-0.07(7) 


-4(3) 


9.1 


2 


1948+3540 


1946+35 


0.717312547254(5) 


7.0592(10) 


54 


340 


6.06(8) 


7.0(10) 


26.0 


9 


1955+5059 


1953- 


[-50 


0.518938338380(8) 


1.3722(15) 


59 


539 


4.0(14) 


6(3) 


21.4 


3 


2002+4050 


2000- 


h40 


0.90506709268(4) 


1.743(6) 


53 


757 


-2.6(10) 


-0.9(62) 


15.3 


6 


2013+3845 


2011- 


h38 


0.230195129550(7) 


8.8514(11) 


50 


602 


-1.52(20) 


-3.7(12) 


15.4 


6 


2018+2839 


2016- 


[-28 


0.557953548298(8) 


0.1478(13) 


42 


282 


-1.08(3) 


-1.6(13) 


30.1 


1 


2022+2854 


2020- 


h28 


0.343402486392(4) 


1.8924(6) 


44 


203 


0.3(3) 


-1.1(8) 


12.8 


14 


2022+5154 


2021- 


h51 


0.529198257306(18) 


3.061(3) 


72 


884 


-0.73(4) 


-4(3) 


9.1 


2 


2048—1616 


2045- 


-16 


1.96157730029(5) 


10.973(9) 


41 


474 


-2.15(20) 


12(9) 


30.1 


1 


2108+4441 


2106- 


[-44 


0.414870558915(12) 


0.095(3) 


52 


679 


-0.02(5) 




9.1 


2 


2113+4644 


2111- 


[-46 


1.01468510534(4) 


0.719(7) 


48 


834 


-0.84(18) 


7(7) 


9.1 


2 


2157+4017 


2154- 


h40 


1.52526635476(9) 


3.434(14) 


61 


1308 


4.0(3) 


8(14) 


9.1 


2 


2219+4754 


2217- 


[-47 


0.538470040747(10) 


2.7634(20) 


57 


800 


0.13(3) 


-1.6(20) 


9.1 


2 


2257+5909 


2255+58 


0.368247257217(8) 


5.7555(13) 


61 


490 


-0.19(5) 


1.7(13) 


9.0 


2 


2313+4253 


2310- 


h42 


0.349433715696(12) 


0.111(3) 


58 


819 


-2.0(6) 


-4(3) 


21.4 


3 


2321+6024 


2319- 


[-60 


2.25648988440(7) 


7.017(13) 


58 


1322 


0.3(6) 


-20(13) 


9.1 


2 


2326+6113 


2324- 


h60 


0.233652031097(5) 


0.3513(8) 


53 


360 


-0.008(10) 


-1.3(8) 


9.0 


2 


2354+6155 


2351- 


h61 


0.94478711093(3) 


16.248(5) 


56 


685 


-0.59(12) 


-16(5) 


9.1 


2 



° 1. Manchester & Taylor, 1981 2. Arzoumanian et al., 1994 3. Backus, Taylor & Damashek, 1982 4. Shemar 
& Lyne, 1996 5. Siegman, Manchester & Durdin, 1993 6. Dewey et al., 1988 7. Manchester et al. 1983 8. 
Newton, Manchester & Cooke, 1981 9. GuUahorn & Rankin, 1978 10. Vivekanand, Mohanty & Salter, 1983 
11. Wang et al., 2000 12. Clifton et al., 1992 13. Manchester et al., 1996 14. Foster, Backer & Wolszczan, 1990 



Table 3. Pre-glitch ephemeris of the Crab pulsar. 



Table 4. Rotation parameters after the 2000 July glitch. 



PSR 


B0531+21 (J0534+2200) 


RAJ 


05'^34™31!972 


DECJ 


22°00'52'.'07 


PEPOCH 


51562.7279 


^ (s-i) 


29.845547780 


V (s-2) 


-3.7457341x10-1° 


i> (s-3) 


1.0161006x10-20 


V (s-4) 


-6.0x10-31 


DM (cm-3 pc) 


56.77 



AP and AP are proportional, also consistent with the ob- 
servations shown in Fig. |^. These agreements strongly sug- 
gest that long-term systematic changes in pulsar periods and 
their derivatives are dominated by pulsar glitches (cf. Lyne 
1996; Johnston & Galloway 1999). 



4.3 The Crab Pulsar 

Observations of the Crab pulsar revealed a glitch during 
2000 July. This glitch was also observed at Jodrell Bank 
Observatory and the pre-glitch ephemeris given in Table ^ 
is base d on Jodrell Bank data. They also determined the 



Epoch (MJD) 


51856.0000 


Data span 


51745.2-51966.6 


^ (s-i) 


29.836059670(2) 


(s-2) 


-3.743460(3) x 10-1° 


u (s-3) 


1.17(2) X 10-20 



Residuals relative to the model in Table ^ are shown in 
Fig. ^ and Fig. ^ shows the observed variations in frequency 
and frequency derivative. We estimate the time constant of 
the glitch, Td, to be 4 days. This value gives phase continuity 
at the adopted glitch epoch and we keep it fixed during the 
analyses. Fitting of the glitch model described in Equation ^ 
using TEMPO gives a glitch size Az^/;/ ~ 24(8) x 10"^, an 
increment in A!>/z> ~ 5(2) x IQ-^, and Q = 0.8, that is 80 
per cent of the jump in frequency decays away on the given 
timescale r^. As may be seen in Fig. ^, there is a persistent 
increase in the magnitude of the slow-down rate; fitting for 
this in TEMPO gives A!>p ~ -48(3) x IQ-i^ s-^. Rotation 
parameters for the post-glitch data obtained using TEMPO 
are given in Table y. 

Glitches in the Crab pulsar are frequent and ha ve rel 
ative sizes in the range 4.7 x 10~^ to 85 x 10~^ (Lyne 



glitch e lpoch to be MJD 51740.8 (2000 July 15). We use Pritchard fc Smith 1993| ; [Wong, Backer fc Lyne 200l| ). Lyne 



these parameters in our analyses. 



et al. (1993) showed that two of the largest Crab glitches 
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-40 



-20 



20 



40 



60 



AP (10 ^ s) 

Figure 6. Difference between observed and predicted period derivative AP, plotted against difference in period, AP. The inset is an 
expanded version of the central region. The data set is the same as in Fig. ^ 



0) 







Period Change (10 s) 

Figure 5. Histogram of pulsar period changes. The Crab pulsar 
and PSR B1737— 30 are not included as they suffered glitches, five 
other pulsars with large uncertainties in AP and AP are excluded 
as well. 



{Au/u ~ 37 and 85 x 10~^) had persistent increases in slow- 
down rate, so tliat tlie Crab pulsar is now rotating more 
slowly than it would have without the glitches. It is clear 



that the 2000 July glitch falls into the class of large glitches 
in the Crab pulsar. 



5 DISCUSSION 

We have developed a pulsar timing system for the 25-m Nan- 
shan radio telescope of Urumqi Astronomical Observatory 
and used it to determine positions and period parameters 
for 74 pulsars. Typical timing residuals over one year are 
few hundred microseconds, giving positions typically to a 
few tenths of an arcsec and pulse periods to 0.01 ns or bet- 
ter. Comparison with earlier measurements has shown that 
long-term period fluctuations are probably dominated by 
recovery from (mostly unseen) glitches. 

A cryogenic receiver system operating in the 18-cm 
band for the UAO Nanshan telescope is currently under con- 
struction. This will enable timing and other measurements 
to be made on pulsars with mean flux density as low as 1 
mjy. Compared to larger radio telescopes, frequent pulsar 
observations are possible with this system, making it ideal 
for projects such as the study of pulsar timing irregulari- 
ties and scintillation. The higher sensitivity of the cryogenic 
system will allow more precise timing and observations of a 
larger sample of pulsars. 
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300 200 
AT=t„-t (d) 

Figure 7. A model for the increments in pulse period and period 
derivative due to a glitch occurring 400 days before the current 
time to- A glitch magnitude At-g = 10~^ and a post-glitch decay 
time of 100 days are assumed. Changes for three different decay 

parameters Q and three different permanent changes in period 
derivative APp are plotted. The groups of three lines of one type 
are for different values of Q and the solid, dashed and dotted 
lines represent different values of APp. The post-glitch values 
are functions of Q and independent of APp, so they overlap for 
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Figure 8. Timing residual of the Crab pulsar showing the glitch 
of 2000 July. 




100 200 
MJD - 51590 



Figure 9. A gliteh of the Crab pulsar at epoch MJD 51741. (a) 
frequency residual Ai/ relative to the pre-gliteh solution and (b) 
the variation of !>. The glitch epoch is indicated by an arrow near 
the bottom of each plot. 
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